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ABSTRACT 

We compute the 2-loop thermal partition function of Yang-Mills theory 
on a small 3-sphere, in the large N limit with weak 't Hooft coupling 
A = gy M N . We include N s scalars and Nf chiral fermions in the adjoint 
representation of the gauge group (S)U(N), with arbitrary Yukawa and 
quartic scalar couplings, assuming only commutator interactions. From 
this computation one can extract information on the perturbative correc- 
tions to the spectrum of the theory, and the correction to its Hagedorn 
temperature. Furthermore, the computation of the 2-loop partition func- 
tion is a necessary step towards determining the order of the deconfine- 
ment phase transition at weak coupling, for which a 3-loop computation 
is needed. 
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1 Introduction 



1.1 Overview and outline 

The interest in the thermodynamics of large N Yang-Mills (YM) theories on a small 
3-sphere S 3 , originates in pQ. In that paper, the AdS/CFT correspondence [21 El E] 
was used to study the thermal behavior of M = 4 supersymmetric Yang-Mills theory 
(SYM) on S 3 , with gauge group (S)U(N). As the correspondence dictates, in the 
large N limit with strong 't Hooft coupling A = gy M N , this behavior can be found 
by studying classical type IIB supergravity. This study showed that the theory has 
a low temperature confining phase, in which the free energy scales as N°, and a high 
temperature deconfining phase, in which the free energy scales as N 2 . The two phases 
were found to be separated by a first order phase transition. 

It is then interesting to ask whether this behavior persists even at small 't Hooft 
coupling. In that limit the conformal M = 4 theory can be studied by using perturba- 
tion theory. However the study needs not be restricted to weakly coupled conformal 
theories, and can be extended to asymptotically free YM theories. As opposed to 
conformal field theories, asymptotically free gauge theories generate an infrared scale 
Aym which leads to the breakdown of perturbation theory in flat space at low tem- 
peratures/energy scales. But on S 3 with radius -Rg3 that bad infrared behavior is cut 
off by the finite size of space, i.e. if Rs-s^ym ^ 1 one can use perturbation theory to 
compute the partition function at all temperatures. 

This has been done in [5], [6], where it was shown that for a wide class of models 
the first order deconfinement transition persists even in the limit of zero 't Hooft 
coupling. The phase transition is signaled by a Hagedorn behavior [7] near the crit- 
ical temperature, where the density of states grows exponentially with energy. This 
remarkable behavior can be ascribed to a zero mode of the gauge field (i.e. a mode 
with vanishing quadratic action) which cannot be integrated out in perturbation the- 
ory. Gauge invariance highly constrains the low energy effective action of this mode, 
resulting in non-trivial self interactions for it even at zero coupling. 

To order A the phase transition is of first order and the critical temperature is 
identical with the Hagedorn temperature. However, in the next order in perturbation 
theory [6], the transition might be of first order and slightly below the Hagedorn 
temperature, or second order and at the Hagedorn temperature. In order to deter- 
mine the order of the phase transition one needs to perform a 3-loop computation in 
perturbation theory. For pure YM this has been done in [S] , where it was found that 
the transition is of first order. For YM theories with adjoint matter fields (including 
M = 4 SYM) the order of the phase transition at weak coupling is still an open 
question and examples which exhibit either behavior can be easily constructed. 

In this work we take a step further towards answering this question, by computing 
in perturbation theory the 2- loop partition function of (S)U(N) YM theory in the 
large N limit with fixed, small 't Hooft coupling on S 3 . We include scalars and chiral 
fermions in the adjoint representation of the gauge group with arbitrary Yukawa 
and quartic scalar couplings, only assuming commutator interactions (though our 
computation can be easily generalized to include non-commutator interactions as 
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well). Our computation follows closely the one of [9] where the same computation 
was done for pure YM. 

As an immediate application we can compute the order A corrections to the Hage- 
dorn temperature and to the spectrum of the theory. We check our computation by 
comparing to the known results for M = 4 SYM, which were computed by spin-chain 
methods in [10J. 

The structure of this paper is as follows. In the next subsection we review some of 
the theoretical background necessary for this work. Section previews our conventions 
for the action, and various mode expansions on S 3 . In section [3] we discuss our 
renormalization scheme, and present the computation of the various counterterms. 
In section [4] we present the computation of the 2-loop diagrams which contribute to 
the partition function. Finally, we present our conclusions in section |5j Additional 
information regarding the computation is collected in the appendices. 



1.2 Theoretical background 

In this subsection we summarize briefly some results related to gauge theories on a 
small S 3 that will be used in this work. A more thorough account of these matters 
can be found in [5], [6]. As stated in the introduction we are interested in calculating 
pertubatively the partition function of large N (S)U(N) YM theory on S 3 with 
additional matter fields in the adjoint representation. This can be done by solving 
the Euclidean path integral on S 3 x S 1 with the inverse temperature (3 as the period 
of S 1 , in the weak coupling regime^} 

In order to use perturbation theory we must fix the gauge. It is convenient to use 
the Coulomb gauge, 

d% = 0, (1.2.1) 

where % runs over the directions of S 3 and the derivative is assumed to be covari- 
ant. This gauge condition fixes spatial gauge transformations but leaves temporal 
transformations unfixed. In order to fix those we impose the following additional 
condition, 

d t a = d t ( -i- / A ) = 0, (1.2.2) 



where we have defined the zero mode of the gauge field a. The additional condition 
implies that a is constant on S 3 x S 1 . The mode a plays a special role in our 
computation since it is the only zero mode of the theorjj^J 



As shown in [5] the Fadeev-Popov determinant conjugate to (1.2.2) can be ab- 
sorbed in the integration measure of a, changing it to integration over U = e i/3a . 
Furthermore, it was also shown that the effective action arising from integrating out 
the massive modes is only a function of U and /3, 

2 For asymptotically free theories this regime is Rg^KyM <C 1- Conformal field theories have a 
dimcnsionlcss coupling so we just set A -C 1 . 

3 We assume here that all scalar fields are conformally coupled so their zero mode is lifted. 
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Z{fi) = J [dU]e~ s '"^ = J [dU] (J [dA^dA^d^e-s^ ^ 1 ^ , (1.2.3) 

where by Aq we mean the part of the zeroth component of the gauge field without its 
zero mode a, S F is the gauge fixed Euclidean action on the sphere, and all the other 
massive fields are collectively denoted as <p l . 

The new integration variable U is related to the Wilson line over the time cycle, 
also known as the Polyakov loop V = ^trPexp(J Q /3 A dt) [TT]. The Polyakov loop 
measures the free energy F q (T) of YM theory in the presence of an external quark 
through the relation (V) = exp(—F q (T)), and is a standard order parameter for 
deconfinement. More precisely, 4tr(C7) is the Wilson line of the zero mode of the 
gauge field around S l . 



The form of the effective action (1.2.3) is highly restricted by gauge invariance. 



To 2 loops order, writing S e ff(U,/3) = S ef l oop + Sl f l r° op , it turns out that 



oo 1 



1.2.4) 



n=l 



Q2—loop 
^eff 



pg 2 



N^fnix) (tr(tT)tr(^") - l) + 



71=1 



+ J2 HU n )tv(U m )tv(U~ n - m )+c.c.-2N) 



n, m=l 



1.2.5) 



where x = e P/ R sz and the function z n (x) is defined in terms of the "1-particle 
partition sums" for N s scalars, 1 vector and Nf chiral fermions on 5* 3 , 



ZnlX) 



ZB 



(x n ) + (-l) n+1 z F (x n ), z B (x) = N s z s (x) + z v (x), z F (x) = N f Zf(x) 

_ f~\ — x+x 2 _ (\ _ 6x 2 -2x 3 _ f\ _ 4a; 3 / 2 (1.2. 



Z,[X 



(1-x) 



3, Zy\X 



(1-x) 



3 J 



Z F \X\ 



6) 



(1-x) 3 



The two and three trace terms in (1.2.5) arise from planar diagrams with 3 index 



loops, while the U independent terms come from non-planar diagrams with a single 
index loopR In the low temperature phase, the three trace terms in the matrix integral 



(1.2.3) contribute only through perturbations around the gaussian model, and can be 



neglected to the order in the coupling that we are interested in. 



The matrix model ( 1.2.4 )-( 1.2.5) can be solved in the large N limit by using the 



saddle point approximation. Below the critical temperature the minimizing eigenvalue 



The form of the non-planar terms is as in (1.2.5) since we assume having only commutator 
interactions, so the corresponding U(l) theory is free. Therefore when setting ti(U n ) = 1 and N = 1 
the 2-loop effective action should vanish. This together with the usual N counting determines the 
non-planar terms. 
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distribution of U is uniform. The matrix integral (1.2.3) can be evaluated at this 
saddle point to give^J 



J U(N) 



-A/3F" 



,\f3fn(x) 



1J i _ Zn ( x ) + \ n pf n (x) ' 



1.2.7) 



where F 2 is the non-planar piece given by, 



x) 



"2 ^ ^ fnm(%)- 



1.2.8) 



n,m=l 



For the SU(N) partition function we divide by the U(l) part (assuming only 
commutator interactions so the U(l) part is decoupled from the SU(N) part of the 
theory) , 



Z, 



SU(N) 



Z U{N) 



n=l 



-z n (x)/n+\f3fn (x) 

z n (x) + \n(3f n (x)' 



1.2.9) 



The goal of our computation is therefore to compute f n (x) and f nrn (x). 

The Hagedorn temperature is defined as the temperature for which the large N 



partition function diverges. This happens when one of the denominators in (1.2.7) 

:e the 1-particle partition sums are mo: 
oo this always occurs first for the n 



or (Jl.2.9 ) vanishes. Since the 1-particle partition sums are monotonic functions with 
and z(l) 



z(0) 



1 denominator in 



the multiplication of (1.2.7) or (1.2.9), so the Hagedorn temperature is given by 
Z\{xh) = 1 • The order A correction to the Hagedorn temperature is therefore, 



Sx H = f3 H X^\, (1.2.10) 

where xh is the Hagedorn temperature of the free theory. 

Other quantities that can be extracted from the partition function are the sums 
of the order A corrections to the spectrum of YM on S 3 x R. These can be obtained 
by expanding the partition function in small x (small temperatures^ 



Zsu { n)= £ * E ^ = 5> f [d{^)+XHx)Y^SE n>i , (1.2.11) 

all states n=4 \ i / 

5 The planar gaussian action vanishes at this saddle point and only contributes to the partition 
function through the integral over the fluctuations around it. That contribution is independent of 
N (up to some irrelevant overall normalization of the partition function) . This is why we have to 
include also the non-planar terms in the partition function around this saddle point, even though in 
the action they are suppressed by a factor of iV 2 compared to the planar terms. 

6 For the U(N) theory there are also a states with energy l/Rs^ and 3/2Rga corresponding to 
the scalar operator tr(0) and fermionic operator tv(ip), respectively. 
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where in the second equality the sum is over states with energy with degeneracy 
d (|), and 5E n ^ is the order A correction to the state % with energy t^t^-. 

By using the state-operator mapping of conformal field theories, the perturbative 
energy corrections of the theory on S* 3 x R are equivalent to the set of 1-loop anomalous 
dimensions of the theory on R 4 . All the theories that we are considering are classically 
conformally invariant and remain so to leading order in A. Therefore we can use the 
partition function on S 3 to compute the sums of anomalous dimensions in flat space. 
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2 Setup of the calculation 



In this section we present some notational conventions and the basic ingredients used 
throughout the computation. 



2.1 The action on S 3 . 

We write the Euclidean action of YM theory on S* 3 x S* 1 with the inverse temperature 
j3 as the period of S 1 . We include scalars $ (a, b, . . . = 1, . . . , N s ) and SO (4) Weyl 
spinors ^ (/, J, . . . = 1, . . . , Nf), all in the adjoint representation of the gauge group 

[/(AO, 

S E = f dt f (KhA-F^F^ + J$ a (- J D 2 + l)$ a + ^ t V J D^ / - 
Jo Js 3 L4 2 

- ^ 2 Q afecd $ a $ 6 $ c $ d + ^g^Y^e^^j] + ^gp^'el**,^]}, (2.1.1) 



ig[A 



* . 



(2.1.2) 



The appropriate veilbein and connection terms for the vector and fermion covari- 
ant derivatives should be understood. We also take the scalar fields to be conformally 
coupled by adding ^<p 2 a to their kinetic terms^j 

Without loss of generallity we take the quartic scalar coupling Q abcd to be cyclic 
in its indices, and the Yukawa coupling p a u to be antisymmetric in J, J. This is the 
most general action that is classically conformally invariant, up to the assumption 
of having only commutator interactions. The reader is referred to appendix [A] for 
further conventions regarding spinors and the restriction of the above action to the 
M = 4 SYM case. 



Next we impose the gauge fixing conditions (1.2.1 ) and (1.2.2). The Fadeev-Popov 



determinant conjugate to (1.2.1) leads to the addition of the ghost action, 



S FP = -J dt J &xto{c%D' 



(2.1.3) 



The resulting gauge fixed action can be found in appendix |B.1| Next we expand all 
the fields in Kaluza-Klein modes on the sphere^} 

A (t, fi) = a a (t)S a (ty, c(t, ft) = c a {t)S a (Q), A(t, ft) = ^ A a (t)V a {Q), 

a a a 

^j(t,n) = J2i>nt)y a (ty, $a(t,n) = x;c(*)<s Q (ft), (2.i.4) 



7 The scalar kinetic term is conformal in 4 dimensions if we add j^TZcj) 2 to the action. We work 
in units for which Rg3 = 1 so the Ricci scalar turns out to be 1Z = 6. 

8 Note that in principle we could have expended the vector field also in dS a (f2), but the Coulomb 
gauge excludes such terms. 
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where S a (Q), V a (Q) and y a (Q) are scalar, vector and Weyl spinor spherical harmon- 
ics. The Greek letter indices collect within them all the relevant angular momenta 
quantum numbers associated with those functions and summation over all of them 
is implied. To write the interaction terms we also need to compute various integrals 
involving three such spherical functions 

£>a/3-y _ I S a £0g1 Quf3-y _ I £ a yP . \75 7 £) a Pl — / "p Q . y/ 3 £7 

Js 3 Js 3 Js 3 

E ap-y = v a • (v^ x v 7 ) , f q/37 = / y^y^s 1 , G aPl = / y a] a^y^v 1 , 

Js 3 Js 3 Js 3 

jjapy = / y a ey ,3 S 1 j a ^7 = / ya\ ^y^ 'S"f K a P~( s = / ^"SPS^S 5 

Js 3 Js 3 Js 3 

(2.1.5) 

Properties of S s spherical harmonics, and the solutions of the above integrals can 
be found in appendix O 



The action in terms of the above expansions and integrals is written in (B.1.4)- 



(B.1.6). Since the fields a a and c a appear in the action only quadratically, it is 
convenient to integrate them out first. This was done in the appendix B.2| and the 
resulting action after integrating out is, 



dtU: { \a\-D\ + (j a + l) 2 )A a + \K{-D\ + ( Ja + 1) 2 )C+ 



+# ta/ (A + e a (j a + \ , (2.1.6) 



2 

»/3 



S 3 = ig dttxL a (j a + \)E a ^A a A fi A 1 - C a ^[A^ 



^]ct>l-iG a ^^ aI [A\^}- 



— 2 F^{p a] ) IJ [f a ^j] + - 2 F a ~ M ^ Ia pUM ]J "} \ , (2-1.7) 



1 r p f 1 - 

S 4 = — g 2 / dttr \ D aiX B Xf3S [A a , 0f] [A\ <p 5 a ] + -B apx B x ^Q abcd (j)a^ 

2 Jo I 2 



naf3\ r)-y5\ f>a/3\ U'ySX 

- h{hJr2) [A^D t AP\[A\D t A 5 \ - jxUx + 2) [€,D t miDM}+ 
+ -j— V M^ ]aI W & j^ J } , (2-1.8) 
where D t = d t — i[a, *]. The other effective vertices which arise from integrating out 



a a and c a do not contribute to our computation as explained in appendix B.2 



9 The integrals H a ^ , J a ^ and K a/3lS needed for the Yukawa and 4-scalar interactions are 
redundant and can be expressed by other integrals: H a ^ — F a ^ J , J a @"i = — F a ^ J , and K a ^ lS — 

QXa0 g\j5 _ 
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2.2 Propagators 

The propagators of the theory can be computed from (2.1.6 
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m{t%M{t))ki) ■ 

We also used the following correlators, 



A«' fcj \t' - t)6 al3 , (2.2.1) 
A^{t'-t)5^5 ab , (2.2.2) 

e^'^'-^rt/j. (2.2.3) 



(D t A%{t') DtAfot)) = 5(f - t)S^S a 5 kj - (j a + l) 2 A^(t' - t)8<#, 
(AC(OtfW)) = -(ti?V)D t <f/g (*))> = AA^'(t' - t)5^ b , (2.2.4) 
(W) A0 6 fc ?(*)) = - t)6°f'6 a 6 kj 6'*- 



-(j a + l) 2 A^(t' -t)5 af3 5 ab . 

The bosonic propagator A(i) is a periodic function with period (3. For < i < (3 
it is defined as, 

\il,kj/.\ _ / e' af ( e-(JA + 1 ) i e (jA + l)t \ \ tf ' fcj ' Icy r> r\ 

where a = a ® 1 — 1 (g> a and a term in the expansion in powers of e lf3a = e l ^ a <g> e _i/3a 
is understood to carry indices by putting them directly on the tensor product ((A ® 
B)u,kj = AuBkj). 

The fermionic propagator 0(i) is anti-periodic with (anti)period (3. We can define 
it over t G (— f3, (3] (and continue periodically) , 



6 A (t) = e« x { (2.2.6) 

l +e -i/3a e /3 e A(JA+2) 



iG(0,/3] 

te (-/3,o] 



where indices should be put on it by the same prescription as in the bosonic case. 

Note that the fermionic propagator, and the derivative of the bosonic propagator, 
are ambiguous at t = 0. We found that this ambiguity only leads to an overall 
normalization of the partition function, which we fix by hand anywaj- 



ii 



2.3 Counterterms action 

As usual in perturbation theory, our 2-loop computation contains divergences. We 
use the regularization scheme of [Bj, which involves cutoff regulators which break 
gauge invariance. This scheme involves multiplying each vector, scalar and fermion 

10 We use i, j, k, . . ., as fundamental representation indices. 
n Say to Z(T= 0) = 1. 
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12 



111 

de- 



line which carries momentum p, with smooth regulating functions R g , s ,f 
the argument of the regulators, M denotes the cutoff regularization scale. We 
mand that all the regulators have Heaviside-like properties: R(0) = 1, R'{0) = 0, 
R(k — > oo) = 0, but other than that they are arbitrary. 

These regulators preserve rotational symmetry on S 3 , but break Lorentz and gauge 
invariance. It was shown by 't Hooft [12J, that gauge invariance can be restored by 
adding the appropriate set of renormalizable local counterterms which may break 
Lorentz and gauge symmetries, while preserving rotational symmetry. To first order 
in A = g 2 N, adding the following counterterms is necessary^ 

Set = A J dt J tx SXJ W j^(Z 0g - Z lg d 2 - Z 2g D 2 )A t + 
+ <$> a {Zf s - Zf s d 2 - Zf s D 2 )<$> b + 



+ MV(Zij<j% + Z*}D t )*j} = (2.3.1) 

A / dttr 8 ^ [A*(Z 0g + Z lg ( 3a + 1) - Z 2g D 2 )A a + 
+ + Zf s j a (j a + 2) - Zf s D 2 t W h + 

+ i^ aI {Z[ J f e a {j a + i) + Z{]D t )r\ 



The notation in (2.3.1) suggests that we take the trace only over the SU(N) 
part of the fields, (e.g. A i ^ = Ai — 4tr(Aj)). Since we assume only commutator 
interactions, the U (1) part of the theory is free and only the SU (N) part of the fields 
contributes to the counterterms. 

For the gauge boson and fermions, imposing gauge invariance is sufficient in order 
for the theory to remain conformally invariant at order A. However, for the scalar 
fields this is not so since their mass terms are gauge invariant and contribute to the 
partition function at order A. As noted in the introduction we want our theories 
to remain conformally invariant to leading order in A. Therefore, we need some 
renormalization condition that keeps the scalars "massless" on S 3 . In practice, we 
set the counterterms Z^ by requiring that the state-operator mapping works (which 
is why we wanted our theories to be conformally invariant in the first place). It 



turns out that the knowledge of the sum of the anomalous dimensions of dimension 

7ab 
J 0s 



2 operators in flat space, is sufficient to determine Z^ completely. 



12 On S 3 we used the prescription which attaches R g ^ s {^f) for vectors and scalars, or Rf y^jf- 
for fermions with total angular momentum j. 

13 Note that we did not write a Zq term to the fermions even though it is allowed. It turned out 
that this term does not contribute in our computation. 
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3 Regularization and counterterms 



In this section we present the computation of the various counterterms in (2.3.1). 
The counterterms Z\ and Z 2 of all fields can be computed in flat space, since from 
dimensional analysis they cannot depend on Rs^. The Zq counterterms on the other 
hand, can depend on the global properties of S 3 . For instance the scalar terms may 
include a dependence on the Ricci scalar, ZQ b (j) a (j) b D Z' 7l(f) a (f)b, and vector terms may 
also depend on the Ricci tensor. 

Following [9j we determine the counterterms by computing 1-loop diagrams in the 
cutoff scheme and demanding that the result match dimensional regularization with 
some convenient subtraction scheme. This guarantees once and for all that the results 
of the 2- loop calculation are gauge invariant. 

3.1 Flat space counterterms 

For the dimensionally regularized theory we work in 3 + d spatial dimensions, and 
extend the coulomb gauge condition to include the d extra components of the gauge 
field. The flat Euclidean action is then, 

5 2 = J d^xtv^AiAi + \djAidjAi + ^AodiAo + d^d lC + 

+ \d^ a d^a + # f V^/}, (3.1.1) 

53 = 9 J rf d+4 xtr|ii i [A i , Aq] - idiA [Ai, A ] - 1^^, A,-]- 

-idiJ[Ai,c] +V t V ) [A ,^j] +V t V[A 4 ,Vi] -i<l>*[AoM- 
-idMAAa] + \^ie{p a] ) IJ [ct>a^j} + l^efiAfa,^']}, (3.1.2) 
S, = ~\g 2 j d d+A xtr{[A , A,] 2 + ^[Ai,Aj] 2 + [A o ,0 a ] 2 + 

+ [A,<p a } 2 + l -q, abcd ct>aMc(t>d}. (3.1.3) 

We define the spinor matrices to satisfy <j^(t u + a u a^ = 2<5 / " / l, with 
fi, v — 0, . . . ,d + 3. We also use the conventions tr(l) = 2, and a = a = 1, which 
are consistent with the d — > limit. The propagators of the fields are, 
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(M", k)M-v, —k)) = Ay(i/, k) = f a ^~*ffi , (3.1-4) 

(A (u,k)A (~u,-k)} = ^ (3.1.5) 

(c(u,k)c\u,k)) = ^, (3.1.6) 

(0>, —k)) = -^-7T 2 5 ab , (3.1.7) 

z/^ + k A 

{ipi{v, k)ip^ J (v, k)) — —^jt^, ■ (3.1.8) 

z/^ + Ar 

In order to compute Z\,Z 2 , we'll compute the 1PI self-energies of the various 
fields. The self energy diagrams that contribute are depicted in figure 1. No quartic 
vertex contributes to these counterterms since such diagrams don't depend on external 
momenta. Also for the gluon self energy the A loop cancels with the ghost loop, so 
we don't have to compute diagrams SEle and SElf. 

Diagrams SEla and SElb were evaluated in [9]. The expressions for the other 
diagrams are (omitting A and group indices throughout.), 



Vector self-energy: — ^(Ai(u,p)Aj(—u,—p)) 



ipi 



r dud 3+d k ky - 2fc j fcJ 

V (2ir)±+ d (y 2 + -u)* + (p-k) 2 ) ' { } 

QF1H-9/V f dud3+dk - ") + * • (P - fcg + 2fc'(fc - pY 

SEld - 2Nj J (l/2 + fc 2 )((w _ I/)2 + (p _ A . ) 2 ) • ( 3 - L1 °) 

Scalar self-energy: — |(0 a (c<j,p)0 6 (— u, — p)) 1FI • 

sab f dvd 3+d k {2p - k) • A • (2p - k) 
(2vr) 4 + d (lo - z/) 2 + (p - k f 

SF2h- ^ ! {2U ~ Vf (1M2) 

^ 2b -~ d J (2^k 2 {{u-vy + {p-ky) ' (3 ' L12) 

SE2c - trf o°V + o b ^o a ) f dv ^ k v{uj - u) + k ■ {p - k) 

SE2c-tr(p p +p p)j {27r)i+d {]y2 + k2){{uJ _ u)2 + {p _ k)2) ■ (3.1.13) 

Fermion self-energy: (ipi(ui,p)ip^ J (uj,p)) 1FI . 



SE2a = -5<* I /0 ,,^ V, ±, , ,y M2 7 , (3.1.11) 



dvd 3+d k 2[{p - k) ■ A}^ 1 - A"[u - v + (p - fc) i a i ] 



SE3a = 25/ / - — ^ — ^ — \ — ^ , (3.1.14) 

7 ./ (27r) 4+d (w - z/) 2 + (p - /c) 2 V 7 

SF o h _ 2 ,J / ^ 3+ ^ q; - y - (p - k) i a i 

SE3h ~ 25 i J P((c-^ + (p-fc) 2 ) ' (3 ' L15) 
SE3c - -2f^V /" ^ 3+d fe ^ ~ v_ + (P ~ A;)^' 
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SE1a SE1b 




SE1c SE1d 




SE2a SE2b 




SE3a SE3b 



i l 

SE3c 

Figure 1: Self energy diagrams contributing to Zl and Z2. Dashed (dotted) lines 
stand for the scalar (A ) field, wiggly line for gluon and solid for the fermion. The 
dashed arrowed line stands for the ghost. 

For the Z 1: Z 2 counterterms we need to extract the coefficients of p 2 ,uj 2 in the 
vector and scalar fields self-energies, while for the fermions those would be the coeffi- 
cients of UjPicr 1 . In order to obtain the cutoff scheme integrals we can just set d — > 
and multiply the propagators by the relevant regulator functions. We will start with 
the dimensional regularization expressions. 

3.1.1 Dimensional regularization in flat space 

The coefficients of external momenta can all be expressed by the following integral, 
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dvd 3+d k 



v 



2mjJ2n—2m—2 



Art 



(27r) 4+d (k 2 + a 2 )(v 2 + k 2 ) n 
1 r(n — m— |)r(m+|) ^ ^ 



rn 



+ 



1 1 , / \ 
- + -ln (7 r) 



-? + l 



(3.1.17) 



Here, p is the regularization scale and e = — d. The subtraction scheme we use 
is to set the square brackets term in 1(a) to zero. The relevant coefficients of p, uo in 
terms of I mn 's , for vector, scalar and fermion self energy diagrams, can be found in 



appendix D.l.l 



Note that l m ,n(0) contains IR divergences, which of course do not exist on S 3 . To 
deal with this we can write 7(0) = 1(a) + [1(0) — 1(a)]- Since the term in square 
brackets contains no UV divergences it must match between the two schemes. So 
we are allowed to use only the IR regulated term 1(a) in order to compare between 
the two regularization schemes. The full expressions for the IR regulated terms after 
applying the subtraction are, 



-1(^(0^)^(0,- 
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(3.1.18) 



In ( + 



(3.1.19) 
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ip i 

„2 



V 
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1 1 



67T 2 



48tt 



167T 2 

f tr (p V + pV) 



In ( ^1 + 



(3.1.20) 
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(3.1.23) 



3.1.2 Cutoff regularization in flat space 



We take d — > in the dimensional regularization expressions (3.1.9 )-(3. 1.16), and 



multiply each diagram by the appropriate regulator functions as described in sec- 



tion |2.3| We now repeat the p and u expansions, including the regulator function 
expansion when necessary, 
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R' (A) (p • k) 
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2M 2 k 2 
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R'(^)p 2 



2Mk 
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(3.1.24) 



After picking up only the relevant coefficients of the external momenta, the in- 
tegration over v is always easy to perform. We then write our results in terms of 



regulator dependent integrals defined in (D.1.21 )-(D.1.23). 



We also include the counterterm diagrams arising from (2.3.1), which contribute 
to the self energies, in the total result (figure 2). 



Figure 2: Counterterm contribution to 1-loop self energy diagrams. 
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(^(w,0)V tJ (a;, 0)^ = 0; 



r (pV f )/ln 



a y 
(3.1.30) 



3.1.3 Flat space counterterms summary 



We solve for the Z's by equating the expressions in (3.1.18 )-(3.1.23 ) with the ones in 



(|3.1.25|)-(|3.1.30|). This determines the flat space counterterms, 
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3.2 Curved space counterterms 



In this section we'll compute the Z counterterms in (2.3.1). As described in the 



beginning of this section the Zq counterterms may depend on the global properties 
of S 3 , so it is necessary to perform a computation on S 3 to determine them. For the 
gauge field counterterm Z 0g we apply our gauge invariant renormalization scheme, i.e., 
we compute a set of correlation functions on S 3 using both dimensional regularization 
and cutoff scheme, and then match the results. We do that by computing the 1-loop 
1PI self-energies of the lowest total angular momentum mode of the gauge field. 

For the scalar field, using this scheme to determine the counterterm Zq s leads to 
the breaking of conformal invariance at order A. We will determine Z^ by demanding 



that our final result is conformally invariant in section 4.3 In this section we will 



only compute the structure of divergences subtracted by Z^, by computing the 1-loop 
self-energies of the lowest total angular momentum mode of the scalar fields in the 
cutoff scheme. 

Throughout this subsection we work in units where \x = -R53 = 1, so the reader 
should not be alarmed when seeing functions which depend on dimensionful quanti- 
ties. These always consistently vanish between the diagrams. As before we start with 
dimensional regularization before moving on to the cutoff regularization method. 



3.2.1 Dimensional regularization on S 3 

To apply dimensional regularization on curved space we let our theory live on S 3 x 
H d+1 with d = — e. We extend the Coulomb gauge condition to involve all components 
of the gauge field other than A , R 



diAi + d a A a = 0. (3.2.1) 

14 In this section A a will denote the R d components of the gauge field, A and Ai mean the same 
as before. The scalars internal indices will be denoted by a,b... 
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Now Ai has a different mode expansion than before since its divergence doesn't 
vanish anymore. The expansions of Ai and A a are, 



A a = Aa a S °> (3- 2 - 2 ) 



aj o >0 



A i= E U a V?+ ■ T^—^ d a A a a diSA . (3.2.3) 
The quadratic part of the action for the gauge field and scalars are then 

Si + S* = J dtd d xtr Qa*(-# - d 2 a + ( Ja + l))A a + \a%-d 2 a + Ja ( Ja + 2)K+ 
+ \^{-d a t - d 2 c+Ja (j a + 2))(5 ab - ja J^ 2) )A a + (3-2.4) 



The additional interaction vertices are spelled out in appendix D.2.1 



For the spinor part a convenient basis for the Dirac matrices of SO(d + 4) ii 
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r^r^r, r a = i®r a , (3.2.5) 



where {r, T a } = and {r a , T b } = 25 ab tQ The T a and T can be taken to be in the 
Dirac representation of SO(d + 2). 

The only terms we will need in the Dirac Lagrangian on R d+1 are then (vielbein 
and spin connection should be understood.), 

C Rd+1 = I tr {iX (f% + f-d a ) A + g\r [A, A] } = 

= [ tr Ux j (r* ® r°r^ + 1 ® r°r a <9 a ) a + #aV ® r°r [A, A] } = 

is 3 

= tr (T^ a + e a (j Q + 1)T) r - igG^T [A\ } , (3.2.6) 
where in the last line we expanded A in S 3 spinor spherical harmonics , 

\ = y a (n)®ip a (y). (3.2.7) 

From the quadratic parts we extract the propagators, 

15 Here a = 0, . . . , d label M, d+1 components. Hopefully there would be no confusion with the 
decomposition of the gauge field in which a labeled d directions in H d+1 , since the other component 
(0) does not participate in the coulomb gauge condition. 

16 We used tr(l) = 1 for the unit matrix in T a space. This is consistent with the d — ► limit, since 
the action should then reduce to the Id action we had after the mode expansions on S 3 . 
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(a a (u, k)a\-v, —k)) = 5 a ?A a (k) = 
(A a {u, k)A^{-v, -k)) = 5 a ^A a (u, k) 



k 2 +j a ( Ja + 2) 
5 a ? 



u 2 + k 2 + ( 3a + iy 



(A:(u,k)A p b (-iy,-k))=8^A a ab (u,k) 



5 aP 



k n k, 



v 2 + k 2 + Ja ( Ja + 2) ^ ab k 2 + j a (j a + 2) 



{<%(u, k)<fi(-v, -k)) = 5^5- ai A a (u, k) 
(if)f(v t k)ip WJ {v, k)) = 5 af3 5j J e a (u, k) 



u 2 + k 2 + (j a + 1) 
5 a(3 5 I J 

k a F" + r\ a 



(3.2.8) 
(3.2.9) 

), (3.2.10) 
(3.2.11) 
(3.2.12) 



Now, besides the diagrams in figure 1 we have additional diagrams which are 
shown in Figure 3. 




SE1g 



SE1h 



SE1i 




f \ 




'OOOOOOOOOOOOOOOOOOOO' '00000000000000000000" ' oooooooooooooooooooo 1 'oooooooofioWoooooooo 1 

SE1j SE1k SE1I SE1m 



■- - 
/ \ 
I I 




SE2d 



SE2e 



SE2f 



Figure 3: Additional self energy diagrams when using dimensional regularization on 
S 3 . Dashed (dotted) lines stand for the scalar (a) field, wiggly (wavy) line for A (A) 
and solid for the fermion. 

We compute the gauge field 1PI self-energy (A a (0)A l3 (0)), with external momenta 
in the directions all set to zero. 

The expressions for diagrams which were not already evaluated in [U] are, 
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/dud d k 
—— jA^QA'fok), (3.2.13) 

" " * tr (9>, fc)9 5 (z/, fc)) , (3.2.14) 
(27r) d+i 

/dud d k 
j—^A s (u,k), (3.2.15) 

To evaluate these diagrams we set the quantum numbers of the external lines to 
(3 = a and the total angular momenta on the external lines to j a = jp = 1. We can 



then sum over all quantum numbers (using appendix C.3) which do not involve the 
propagators. This usually leaves us with simpler integrals due to the selection rules 
imposed by the 3-j symbols. After solving the integrals the summation over the rest 
of the quantum numbers will be regulated by d. 

For example, in SElc the selection rules constrains the summation to j 7 = jg 
and jry > 1. The angular momentum sum over all quantum numbers but j 1 gives: 
^(jtaBcps* = 2R 2 c (j 1} 1, j 7 ) = MMWrtg) The integral is easily performed and we 
get, 



^yfU £ + 2)(; 7 + ly- 1 = ^^S(c(-i - v - c(i - d)) 

vr ( 4?r ) a tt ( 47r ) 2 



J7 = 1 

= -^([i + ^W-| + l] + 7 + ^) + 0( £ ). 

(3.2.16) 



Other diagrams can be evaluated in a similar way. The results diagram by diagram 



are given in appendix D.2.1 After employing our subtraction scheme (see below 



(3.1.17)), the total result for the 1-loop 1PI self-energy of the gauge field is, 



Q-N s - AN f 9 + N s + 16N f 10 
°DR = ~ 2 h + ^ - ^C(3). (3.2.17) 

3.2.2 Cutoff regularization on S" 3 

For the gauge field we obtain the cutoff regularization expressions by taking the d — > 
limit in the dimensional regularization expressions, and multiplying with the appro- 
priate regulator functions. We also compute the 1PI self-energy diagrams (figures 
1 and 3) of the scalar fields, (<j)^(0)m(0)) with (3 = a and j a = jp = 0. The ex- 
pressions for those (using the notation from the last section for the propagators, and 
suppressing the regulators) are^ 



17 Qab(cd) — Qa/icrf _|_ 



21 



SE2a - gS ab C a ^ 5 C^ S I — A 7 (z/,0)A d >,0), 

2n 

SE2b = 25 ab B a ^B^ s I — z/ 2 A 7 (0)AV,0), 

2n 



SE2c 
SE2d 
SE2e 
SE2f 



-tr (p aj! p b + p bj! p a )F^ a F s ^ [ ^tr(e 5 (z/,0)6 7 (-z/,0)) 

J 2n 

r\{ab)cc . _ _ __ . r j 

v_ ^ B cfis B &fi + B ^s B s^ / ji-A 7 (z/, 0), 

/J 
— A 7 (0), 
27T 

_ 2S ab B af36 D ^5 j _Za 7 (z/,0). 

27T 



(3.2.18) 
(3.2.19) 
(3.2.20) 
(3.2.21) 
(3.2.22) 
(3.2.23) 



Note that SE2e has a 5(0) divergence. This term exactly cancels with a similar 
term in SE2b. We express the resulting angular momentum summations with the reg- 



ulator dependent integrals defined in (D.1.21 )-(D.1.23), by using the Euler-Mclaurin 



formula 18 For instance for SElc (3.2.16) we have, 
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(3.2.24) 



The rest of the expressions are given in appendix D.2.2 We denote the total result 
for the self-energies by Geo and Sco for the vector and scalar respectively, including 
the counterterm diagram contribution in the total answer (figure 2), 



Geo = -12 (Z 0g + AZ lg ) + 4M 2 (2C 200 - 6C 2 100 + N s (C 2 020 - 3C 2 010 ) + 4iV / C 2 002 ) + 



(6-N s - 4iV / )(ln(M) + 7) + In 



A 8 



A 2 A Ns A 4Nf 



) + 



+ 



71 + N s - 5N 



f 



6 



- 10C(3) 



(3.2.25) 



Sco = -2^ t + M 2 [(Q^ cc - 25 ab )C% 10 - A5 ab C\ m ] + 
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(ab)cc 



48vr 2 



tr(p at p + p V 
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47T 2 M 2 C 2 002 + - - - (ln(4^oo 2 M) + 7) 



(3.2.26) 



18 For /(")(oo) = the formula is £~ /(*) = Jo°° + §/(°) " tV'(°) + 75q/ (3) (0) + ■ ■ • . 

We never need more terms in the expansion since each derivative lowers the degree of divergence by 
one. 
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3.2.3 Curved space counterterms summary 

We extract the Z 0g counterterm by comparing (3.2.25), with (3.2.17), 



7 9 



-C 2 200 - 2Cl 00 + — C 2 020 - iV s C 010 
3 2 2 3 2 



5 N s 

187T 2 457T 2 J 



M 2 C 9 002 



1 

360tt 2 



7 . i Aioo 
+ „ In 



30^ 



A 



200 



15 



2 
5* 



2 



(2^| + N s r 2 ) + 



(3.2.27) 



The divergent piece of is determined by (3.2.26). The finite piece on the other 



hand is arbitrary if we only demand gauge invariance and we parameterize it as, 



7 ah 
J 0s 
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—JC^ + C 2 tr (p<V + pV) + C 3 Q (af,)cc ). 



(3.2.28) 



The constants Ci, C 2 and C3 will be determined in section 4.3| by demanding 
conformal invariance to order A. 
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4 Perturbative calculation of the effective action 

In this section we present the computation of the large N effective action S e ff(U) 



defined in (1.2.3), to 2-loops in perturbation theory. 



4.1 2-loop vacuum diagrams 

The 2-loop diagrams which we have to evaluate are depicted in figure 4. 




L5b 



Figure 4: The 2-loop diagrams obtained after integrating c and a. Vertices with a 
little circle correspond to terms we got after integrating the linear a vertex. 

Diagrams Li, L 5a and L 5 b were already evaluated in [9j. The expressions for the 
other diagrams that contribute in the planar limit are, 



L 2 = -2f3g*N s C a ^C^ / ^A^(t)A{ fe (t)A^(t), 

Jo 

U = l(3g 2 N f G a ^G^ [ dt{e^(-t)ef(t)A^(t)+ 
* Jo 

+e2(-f)ejy(t)A*(f) 



(4.1.1) 



(4.1.2) 
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J6a - 4 



l -f3ghr{p^p a )F^F a ^ J dt [e»(t)Af(t)ef(t) + e^(t)A^(t)ef(t)] , 

(4.1.3) 

^Pg 2 Q aabb (b^B^ + B afh <B^ A2'(0)A^(0), (4.1.4) 



Lab = tfN. . (D T A«(0)D r Af (0) + & + l) a A«(0)AJ*(0)) , (4.1.5) 

L 7 = f3g 2 N s D a ^B^A^(0)Af(0), (4.1.6) 



pa/37 P/3a7 / 

^ 2 % 7 (j 7 + 2) ( e «(°) e ?(°) + e «( ) e ?(°)J • 



(4.1.7) 



Each diagram has three index loops, and the short-hand indexing on the propa- 
gators indicates how a term in their expansion in powers of e ifia ® e~ il3a should be 
grouped into traces. The notation A u means that the left side of the tensor product 
in the expansion belongs to the index loop i and the right side to the index loop j. 
Also, the propagators above should be understood to be multiplied by the appropriate 
regulators. 

We can sum analytically over all angular momentum quantum numbers, except the 
total angular momentum, by using the identities derived in appendix |C.3| We then 



have to expand the propagators in powers of e l ^ a (g> e i/3 ° and extract the coefficients 



f n {x) and f n m(x) defined in (1.2.5). The results are expressed as sums over the total 
angular momentum numbers constrained by some selection rules from the integrals 
over 3 spherical harmonics. 

In all diagrams we find that we can write f n (x) = f\ i+ (x n ) + (— l) rt+1 /i i _(x n ). 
The expressions for those are long and not very illuminating. They are summarized 
in appendix [Ej 

4.2 Counterterm vacuum diagrams 

The counterterm diagrams which contribute to our computation are depicted in figure 
5. For the gauge field and scalar we find that fni x ) — fi{ xU ) = f ct {x n ). After 
summation over angular momentum we obtain 
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ff{x) = 2 (Z 0g f(x) + (Z lg - Z 2g ) gi {x)) , (4.2.1) 
ff(x) = (k(x)Z™ + 9l {x)Z? s - ( 9l (x) + k(x)) Z£) , (4.2.2) 

where, 



k(x) = f(x) = k(x) + ln(l-x), 9l (x) = 6k(x) 2 . (4.2.3) 



19 In deriving these expressions we have used the equation of motion of the propagator to determine 
the coefficient of Z 2 , e.g. in the bosonic case (—D^ + (j + l) 2 )A(t) = S(t). In principle we should 
have added the 5(0) term but this term cancels with the corresponding diagram which involves the 
U(l) part of the fields. 
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For the fermion we find fn(x) 
momentum we obtain, 



-l) n 1 f\{x n ). After summation over angular 
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y ^ ^ ✓ 
Lfcf L2cf L3cf 

Figure 5: Counterterm vacuum diagrams which contribute to the 2- loop computation. 

Plugging in the expressions for the counterterm coefficients we have computed 



(3.1.31)-(3.1.36), (3.2.27)-(3.2.28) we obtain, 
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(4.2.7) 

In principle we are done and all we have to do is to plug the expressions for f n 



and f nm we found (E.l.l 



-(E.1.13), (E.2.3)-(E.2.12) summed with the counterterm 



expressions, (4.2.5 )-(4.2.7) inside the partition function formula (1.2.9) and perform 



the angular momentum summations numerically. 

However our regularization involved several subtleties, and it is useful to check 
explicitly that all regulator dependence cancels between the counterterm vacuum 
diagrams and the 2-loop vacuum diagrams. 

In order to see this cancelation we first extract the regulator dependent part 
of the 2-loop diagrams by expanding asymptotically in large angular momentum 
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as was done in [9]. The expressions we found for the regulator dependent part of 
each 2-loop diagram are listed in appendix |E.3 One can readily check that indeed 



ft + ft + Yli f r if+ an d ff + Yli flT- are re g ma tor independent. This is a useful 
intermediate self consistency check on our computation. 



4.3 Numerical evaluation and checks 

As explained in the introduction there are two immediate applications to our calcu- 
lation. One involves computing the order A correction to the Hagedorn temperature 



(|1.2.10|) of the theory. The second involves computing the order A corrections to the 

4 SYM case in pi)] by using 



spectrum of the theory on S 3 , (1.2.11). 



Both have already been computed for the M 

11 tern 

24A ln(x) 



spin-chain methods. Their result for the small temperature expansion 'n 2 ^_ 
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while the correction to the Hagedorn temperature is, 
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A ln(x H ) 
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(4.3.1) 



(4.3.2) 



Using our computation method we can perform this expansion and compute the 
perturbative correction to the Hagedorn temperature numerically. As a self con- 
sistency check we indeed found that the numerics do not depend on the choice of 
regulating function. 



However before we compare our computation to (4.3.1) and (4.3.2), we must fix 
the arbitrary constant piece in (4.2.6). Recall that the source of the constants Ci, C2 



and C3 in ff(x) is the mass renormalization of the scalar fields (3.2.28). Since a mass 



term for the scalars is gauge invariant, these constants are arbitrary. In particular to 
order x 2 our computation yields the following partition function, 



Z(x) 



X< X 1} + (N. + 1) (JV.(48Ci - 33) + (96C 2 + ll)tr(pV) + 



2 ' 48vr 2 
+ (96C 3 + l)Q aabb ) + 6Q aabb ) 



x 2 + 



(4.3.3) 



20 Thc JV = 4 partition function in (4.3.1 1 is adopted to our conventions. Our coupling constant 
is defined to be twice the one in [TU] . 



27 



If our theory is conformally invariant to order A, then the coefficient of x n ln(x) in 



(4.3.3 ) is the sum of the anomalous dimensions of operators with classical dimension n 



in the flat space theory (see below (1.2.11)). In particular, the operators of dimension 
2 are, 



Of = tr(0 a 6 ). (4.3.4) 
We computed the anomalous dimension matrix of those operators in the flat space 



theory corresponding to (2.1.1 
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[3N S (N S + 1) - (N s + l)tr(pV) - Q aabb ] 



(4.3.5) 



dim =2 



Comparing the coefficient of the x 2 lia(x) in (4.3.3) to (4.3.5) we can solve for C±, 
C 2 and C 3 , 
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96' 
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(4.3.6) 



We evaluated the angular momentum sums in the 2-loop diagrams numerically up 
to order x 6 (we could easily have reached higher orders with more computer time but 
it was not needed in this work), 



21r 



The assumption of having only commutator interactions is implicit in (4.3.51 since only the 
combination Q aabb appears in it. Under this assumption the two possible combinations are related 
by Q afcQb = -iQ aabb . 
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(4.3.9) 



One can plug those in (1.2.9) and see that for M = 4 (see appendix |A|) we repro- 



duce (4.3.1) and (4.3.2). This is a very non-trivial check of our computation. 

It is also important to have a check for the generic theory which is not conformally 
invariant, in order to make sure that conformal invar iance is preserved to order A by 
our renormalization scheme. To do so we computed the anomalous dimensions of 
dimension 5/2 and 3 operators, 
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in agreement with the results derived from (4.3.7)-(4.3.9). This gives a check that 
by applying our renormalization scheme for the scalar field mass terms, the general 
theory is indeed conformally invariant to order A. 

In the pure YM and Af = 4 SYM theories ATh = tt 2 8Th/Th was found to be 
positive and rational. In the more general cases that we considered, we found that this 
quantity can be either positive or negative and that it is not necessarily rational. For 
example, consider gauge theories with N s +Nf = 10, Q aabb = N s (l—N s ) corresponding 
to tr([(f> a ,<j) b ] 2 ) self interaction, and tr(p a p a t) = N s Nf Yukawa couplings. Strangely, it 
turns out that all these theories have ^(A = 0) = 7 — 4y / 3. The results for ATh for 
these theories are, 
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For pure scalar theories with the same self interactions as above, ATh becomes 
negative for N s > 2 but is not necessarily rational. For instance, 
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For generic couplings the results are of course irrational, but as we see in the ex- 
amples above, rational results are not so difficult to construct. It would be interesting 
to understand the reason for ATh to be rational. 
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5 Conclusions 



In this work we have computed the full large N partition function of SU(N) YM 
theory on S 3 in perturbation theory, to 2-loops order. We included scalar and spinor 
fields in the adjoint representation of the gauge group, with arbitrary commutator 
interactions. The computation involved expanding those fields in Kaluza-Klein modes 
on S 3 . Then, instead of solving momentum integrals, we had to perform sums over the 
angular momentum of those modes and derive various identities in order to perform 
these sums. Furthermore, we used a cutoff regularization scheme which is not gauge 
invariant. To deal with this we added non-gauge invariant counterterms to the theory 
and demanded that their value precisely compensates for the non-gauge invariance 
introduced by the regulator. 

Furthermore, we demanded that our theory is conformally invariant to order A, 
so that we can relate the perturbative corrections to the energy spectrum on S 3 x R, 
to the sums of anomalous dimensions on R 4 , by using the state-operator mapping. 
Though generically our theories have a non-zero beta function, it contributes to the 
partition function only at higher orders. The only obstruction to achieving confor- 
mal invariance at order A is the generation of mass terms for the scalar fields. By 
demanding that the state-operator mapping works for dimension 2 operators, we can 
fix the mass renormalization of the scalar fields such that the theory is still conformal 
at order A. 

Our computation agrees with known results for the 1-loop anomalous dimensions 
of M = 4 SYM which were computed in [10J. This is a highly non-trivial check on 
our computation which was done by a very different method. Our results can be used 
to numerically compute the partition function in the confinement phase of the large 
N gauge theories which we considered. Furthermore, the order A correction to the 
Hagedorn temperature can be computed very easily for those theories. 

This computation also has a direct application towards determining the order of 
the deconfinement phase transition. As shown in [5] this involves a 3-loop compu- 
tation, and the 2-loop computation of fu(x). The latter can be computed from our 
results for the f^(x) 7 s gjgt-pjjgt. 
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A Additional notations 



Having enough indices as it is in our computation, we use matrix notation for spinors. 
Our notation for the Weyl spinor matrices is such that, 

^ = (l,ir*) , ^ = (1,-^), (A.0.13) 

where r l are the usual Pauli matrices. Also we use the charge conjugation matrix in 
the Yukawa terms, 

e = ir 2 . (A.0.14) 

In order to apply our computation to M = 4 SYM, take N s = 6 and Nf = 4. 
Furthermore, in that case the Yukawa matrices satisfy, 

pV + PV = 2rt 4X 4, (A.0.15) 
while the quartic scalar coupling satisfies, 

Q abcd = 25 ac 5 bd - 6 ad 5 cb - 5 ab 5 dc . (A.0.16) 

B The gauge fixed Yang-Mills action on s 3 
B.l The action 



Imposing the gauge conditions (1.2.1), (1.2.2) and adding the conformal coupling for 
the scalar fields, the Euclidean action on S 6 can be written as Se = S g + S s + S/ 
where, 



s. 



J dt J tr l-^A(D 2 t + d 2 )A* - ^A d 2 A - Jd 2 c+ 
+ igD t A l [Ai, Aq] - ig[A\ A ]d t A - igd^A', A*]+ (B.l.l) 

where we define D t = d t — [«,*]• For the pure conformally coupled scalar (kinetic 
plus quartic interaction), 

- ig[A , $„] A$ a - ig[Ai, $ a ]^$ a (B.1.2) 
-\g\A h $ a ] 2 - l -g 2 [A , $ a ] 2 - ^g 2 Q abcd ^ b ^ d 
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and for the fermions (kinetic plus Yukawa interaction) we write, 



S f = [ dt [ tx{i¥ I (D t + a i d i )^ I + 
Jo Js 3 

+ 0*t / [4o,*j]+0*M4,*i]+ 



(B.1.3) 



In the above the spatial derivative is assumed to be the covariant one on S 3 
including the appropriate connection terms for vectors and spinors. 

Using the properties listed in appendix O we can write the action in terms of the 



expansions of the fields in Kaluza-Klein modes (2.1), and in terms of the spherical 



harmonics integrals (2.1.5) as = S2 + S3 + S4, with 
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3xUx + 2) 



+ (d^d^ sx - D aSX D^ x \ A a A l3 A~<A s + B apx B x ^ s [a a , 0f][a 7 , <p s a ] + 



+D a ^ x B Xl3S [A a ,^}[A\(f) S a } + -B af5X B x ^q, 



(B.1.6) 



B.2 Integrating out a and c 



The fields a and c appear only quadratically in the gauge fixed action (B.1.4)-(B.1.6) 
and it is convenient to integrate them out first. Their propagators are, 



W)<#(t)> 
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(B.2.1) 
(B.2.2) 
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From loops of a's and c's we get the effective vertices, 
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The other diagrams arise from strings of a's using the linear a vertices in (B.1.5), 
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^jjafSXj^ySX r/3 
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2,2 



h{jx + 2) 

j^a.pXp~f&X r(3 

jx(jx + 2) J 



dttr ([A a , D T AP]{ij;j, V> t7/ }) , 
dttr ([<!>:, D^M,^ 1 }). 



(B.2.5) 
(B.2.6) 
(B.2.7) 
(B.2.8) 
(B.2.9) 
(B.2.10) 



The 5(0) divergence in (B.2.3) and (B.2.4) is an artifact of the Coulomb gauge 
and arises since a and c are not dynamical fields. We see that those vertices cancel by 
adding to them the expressions obtained by contracting the two covariant derivatives, 
in (B.2.5) and (B.2.6) (see (2.2.4)). Therefore we can forget about 5(0) factors for 
diagrams which involve those vertices. 

Also the vacuum diagrams which arise from one insertion of any of the vertices 
(B.2.8), (B.2.9) or (B.2.10) vanish for kinematical reasons. After the self contractions 
the summation over spherical harmonics (see (C.3.4), (C.3.7) ,(C.3.8)) vanishes unless 
jx = 0, which can not be since this is the angular momenta of a A , (Ao) whose zero 
mode has been factored out. 



C Properties of S* 3 spherical harmonics 
C.l S 3 spherical harmonics 

In this section we list some properties of spherical harmonics that we used. Due to 
lack of space this is far from a complete review, and the interested reader is referred to 
[S], [32] for scalar and vector spherical harmonics, and to [14J,[T5] for spinor spherical 
harmonics. 
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Spherical harmonics sit in representations of the isometry group of S 3 , which is 
5*0(4) ~ SU(2) x SU(2). The scalar spherical harmonics transform in the (j/2, j/2) 
representation, where j > and —j/2 <m,m'< j/2. 

They form a complete orthonormal set of eigenfunctions of the corresponding 
laplacian on S 3 , 

and satisfy the conjugation relation, 



(S™' n )* = ^-l)m+n S -m,-n_ (C.I. 2) 



The vector spherical harmonics sit in the ((j + e)/2, (j — e)/2) representation, 
where e = ±1 and j > 1. This set of functions satisfies the eigenvalue equations, 



v 2 vp n = -(j + iyvj; 11 , 

VxV™'" = -e(j + l)V™/\ (C.1.3) 

V-V™J n = 0, 

and the conjugation relation, 



The spinor spherical harmonics are two-component Weyl spinors, and transform 
in the representations ( J '~^ £ ^ 2 , J ~( 1 + e )/ 2 ) for e = ±1 and j > 1. They can be found 
from the scalar functions by tensor multiplying a (j/2, j/2) representation, with a 
(1/2,0), (or (0,1/2)) basis spinor with the right Clebsch-Gordan coefficients. The 
properly normalized result is 



(C.1.4) 



, ; m,n _ J_ I y / j/2 + mS™ 1 1/2 ' n 



* \ JJj2^nS™X 



y:>- _L_ I -0/2 -m + 1/25, 



m— 1/2, n 



(C.1.5) 



^ I v / J 72 + m + l/2 l S. m+1/2 ' n 



From the properties of the scalar functions one can check that those are normalized 
correctly 

/ y™ e M y™f = Sjji5 m , m >5 njn iS £je >, (c.i.6) 

«/ s 3 

and that under conjugation those obey 

(y£ n Y = (-l)m+n-e/2 £ (^T") . (C.1.7) 
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The spinor functions satisfy the eigenvalues equatio 



a d^; = e(j + -)y 3 ; . 



(C.l. 



In the text we collect all the angular momentum indices of a scalar, vector 
or spinor spherical harmonic by a single Greek index. A barred Greek index on 
a scalar or vector spherical harmonic means the corresponding function is conju- 
gated, while by a barred Greek index on a spinor spherical harmonic we mean 

ya — ^_^m a +n a -e a /2 /y-m a ,-n a ^ 

C.2 Spherical harmonics integrals 

While expanding our fields in Kaluza-Klein modes on S* 3 we defined integrals over 



three spherical harmonics (2.1.5). These integrals can be expressed as a reduced 



matrix element times the appropriate 3-j symbols, 
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(c.2.1; 
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QOLp^f 
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2 
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(j ol -(l-e a )/2) (jp-{l-e )/2) j 1 + 



2 2 2 

n a np n 7 

(j a -{l+e a )/2) (jp-(l+ep)/2) j 7 - 

2 2 2 

n a n^g n 



7 



(C.2.2) 



This can be related to the fermionic integrals listed in (2.1.5), by using the con- 



jugation relations of the spinor spherical functions (C.l. 7). 



We computed the reduced matrix elements which were not already evaluated in 
[S] explicitly, 

22 The inclusion of the appropriate vcilbein and spin connection should be understood. Also recall 
that in our convention for spinor matrices the a % are anti-hermitian, hence the real eigenvalue in the 
above equation. 
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R H +(x,y,z) 

R G 2( x >y, z ) 




(z + l)(a - z)(a + 1) 



(z + l)(cr - z)(cr - y) 



(a - -y)o-(o-+ 1) 
(a - a;) (or - y)(& - z- l)(cr - z) 

(*+l) 

o--z)(o- + l)(o--3/)(o--j/ + l) 



+ 

where a = (x + y + z)/2 and <r = (x + y + z + l)/2 are integers and, 
= Rh++ = Rh — , Rh- = = —Rh-+, 



R-go — R&+++ ~ R& — ' R&i — R& — 



R&2 — Rg+- 



R. 



G-+-1 ±L G2 



~R(3++-i 

R G 2( X v) = - r g-++ = - R - 



(C.2.3) 



(C.2.4) 



C.3 Performing angular momentum sums 

We computed various sums over spherical harmonics quantum numbers, which appear 
in the expressions for our diagrams. We list those identities which were not already 
derived in [HF^l 



Y^g^g^ = (-l)^R^^V a J p JjR^^fo,j a ,j y ), (C.3.1) 

m's 

Y^FtfUpcfo = (-iy^^ +e -^Rl ea£p (j a ,j p ,j 7 ), (C.3.2) 

m's 

23 Translating the notation in [3] to ours: i?2 — * Rc > ^3e Q e^ ~~ ► RDe a ep, and i?4 (t , E(j e 7 — > REe a e3e~- 
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m's 
n's,e' s 

m's 
m's 



(C.3.3) 
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D 1-loop self energy diagrams 

D.l Flat space self-energy diagrams 
D.l.l Dimensional regularization 

We write the flat space 1PI self-energy diagrams (figure 1) which were computed in 



dimensional regularization, in terms of the integrals I m>n = I m , n (0) (3.1.17) 
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D.1.2 Cutoff regularization 

We used the following definitions of regulator dependent quantities in the main text, 



A nkl M\ _ r . qR^R^R^q) 



In = / dq 1 gv " ^ (D.1.21) 

i r°° 

B s { 9 = — / dqR' f (q)R s>g (q), (D.1.22) 

1 /"°° 
«/ o 

1 = — jf dqqR n g (q)R k s (q)R l M- (D.1.24) 
The cutoff scheme results for the flat space self energy diagrams in terms of the 



regulator dependent functions (D.1.21)-(D.1.24) are 
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Scalar p 2 expansion: — ^(0 a (O,p)0 6 (O, — p)) p f 7 - 
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D.2 Curved space self-energy diagrams 
D.2.1 Dimensional regularization 

When applying dimensional regularization on S 3 we get additional vertices included 



in the action (B.1.5), (B.1.6), coming from the fact that the gauge field has additional 



components. We needed the following additional cubic vertices, 
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- . p M0 y , A^)d t d a A a a - . y d a Al\d t A a + D°f»[dJ&, A a }Al+ 
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J7IJ7 T" 2) 
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-A^fc^ - j Q ( Ja + 2)5^[C, d a Al]4>l- 



jfi(jfi + 2)j J (j 1 + 2)' 
-B^iAt^Wl). (D.2.1) 



The additional quartic vertices which were necessary for the computation are, 

% = - l -ghi (5^ A [^V^[AVl+^^1^^[^^Hl+ 

+ C^C^[A a ,d a A^}[d b AlA s ] + — [A a ,^[^,^4g]+ 

J7U7 + 2 >JS{J8 + 2 J 

+ S^ x 5^ A [5 a ^,0S][^,^]) . (D.2.2) 
Above we used an additional spherical harmonic integral, 



B a01 = — — / dS a ■ dS S^. (D.2.3) 

The solutions to the vector self-energy diagrams are, 
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Nil 'v 1 
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7 -|^+0(c), (D.2.5) 



SElk=^ + 0(e). (D.2.6) 

The other diagrams in figures 1 and 3 were computed in [9] (see equation (3.59) 
in that paper). 

D.2.2 Cutoff regularization 

These are the expressions for the self-energy diagrams on S 3 in the cutoff scheme 
diagram by diagram, 
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The other diagrams in figures 1 and 3 were computed in [9] (see equation (4.14) 
in that paper). 



E 2-loop vacuum diagrams 
E.l Planar diagrams 



By expanding the propagators (2.2.5), (3.1.8) in the diagrams (4.1.1 )-(4. 1.7), and 



using the identities in appendix C.3, we extract the f n (x) terms in the effective action 
defined at fll^.Slf 11 
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24 In this section we use a, 6, . . . = 1, . . . , N s for the scalars internal indices. 
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E.2 Non-planar diagrams 

Define the functions, 

}C nm (y : z;w) = [y n + z n } [w m - (yz) m ] 

C nm (y, z, w) = (yz) n (y m + z m ) + (zw) n (z m + w m ) + (wy) n (w m + y m ) . 
In terms of those the fnm s are > 
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The summations over n, m can be performed explicitly, and in particular one has, 
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In terms of those the expressions for F 2 defined in (1.2.8) are 
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In the above, when they were not written, the arguments of the reduced matrix 
elements are R# = R#(a,b,c). 



E.3 Regulator dependent part of 2-loop diagrams 



The regulator dependent parts of the diagrams are obtained from the /j s in (E.l.l )- 



(E.1.13) by expanding the summands in large a and picking up the terms which 
diverge when we take the cutoff to infinit}|^| The results in terms of the functions of 
temperature defined in (4.2.3) are, 



reg 
LI 



j oo 



f{x)a + 



S gi {x)-7f{x)l 



+ 



a=l 

gi (x)-4f(x) 
10M 



R 



a 
a 



+ 



(E.3.1; 



JL2 — o_9 2—t 



3.^1 2 ^(s) + » (ll « fl '© 8 '(E l + 

a=l 



gi (x)-4f(x) 
20M 



^ ' If) K (If) + If ^ s (If) ^ (If 



(E.3.2) 



2 , a-l\ + 2 9l (x)-3f(x) 



0=1 

" ~ / f/ ~ 2 A j^ii ( a ~ 2 



5M 



i2/ 



Mir 



M 



+ 



M 

N f 2x1(1+ a;) v^/l / a \ 



7T 2 (1 -x) 4 i — ' I a 

v ' a=l 



M 



1 :R ' " 



2M 9 \MJ } \ M 



R' f 



(E.3.3) 

1 



rreg 
JL4+ 



tr (p s V a ) 

" 4^ 



a=l 
1 



f \ M 



9i{x) 
3M 



(E.3.4) 



+-R f ( a -^)R"A a —^ 
M \ M I f \M 



1 1 M ) " J V M 

(E.3.5) 



+ 



25 Actually the RR' terms are regulator independent since / R{q)R'{q)dq= - \ but we included 
those anyway. 
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Or in integral form using the Euler-Maclaurin formula and the regulator dependent 



integrals (D.1.21 )-(D.1.24) 
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